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Abstract 

Using a positive semidefinite operator technique one deduces exact ground states for a zig-zag 
hexagon chain described by a non-integrable Hubbard model with on-site repulsion. Flat bands 
are not present in the bare band structure, and the operators b}i,(t introducing the electrons into 
the ground state, are all extended operators and confined in the quasi ID chain structure of the 
system. Consequently, increasing the number of carriers, the -6^,0- operators become connected 
i.e. touch each other on several lattice sites. Hence the spin projection of the carriers becomes 
correlated in order to minimize the ground state energy by reducing as much as possible the double 
occupancy leading to a ferromagnetic ground state. This result demonstrates in exact terms in a 
many-body frame that the conjecture made at two-particle level by G. Brocks et al. [Phys. Rev. Lett. 
93,146405, (2004)] that the Coulomb interaction is expected to stabilize correlated magnetic ground 
states in acenes is clearly viable, and opens new directions in the search for routes in obtaining 
organic ferromagnetism. Due to the itinerant nature of the obtained ferromagnetic ground state, 
the systems under discussion may have also direct application possibilities in spintronics. 
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I. INTRODUCTION 



A. Conducting chains 

Conducting chains have been extensively studied in the past period, due to several rea- 
sons. First, periodic chains holding different type of cells are genuine objects for applications 
in device design techniques at the level of molecular nanotechnologyi. As such, these sys- 
tems represent a starting point in creating functional devices, materials and components on a 
1-100 nanometer length scale which can lead to new routes in realizing functions of practical 
interest like field effect transistors, electroluminiscent diodes, nanocatalysis, hydrogen stor- 
age, etc.-. Second, several type of periodic conducting chains are in fact organic materials. 
Their nature, besides the fact that potentially allows developments that can lead to con- 
densates on a plastic (as plastic ferromagnetism for example^.), has exceptional qualities for 
new developments in electronics, for example in the direction of organic optoelectronic and 
field effect transistor components^, high performance transistors and circuits made based on 
soluble materials^, biodegradable devices used for controUed-release drug delivery inside of 
the human body^, or electronic components on plastic providing flexible electronics^. Third, 
conducting chains exhibit quite interesting properties and phases important in advanced 
technological applications (as spintronics) , for example insulating, conducting, half metallic, 
paramagnetic or ferromagnetic behavior-. Furthermore, these systems have the virtue of 
having fundamentally different ground states which can be tuned by external parameters 
like external fields or site selective gate potentials^, opening new routes for the design of 
valves, switches, or control devices. Finally, the majority of conducting chains being de- 
scribed by non-integrable models, represent genuine challenges to theory and necessitate the 
development of rigorous techniques to describe the physical properties of these systems. 

The chains under consideration are built up in fact from periodic arrangements of rings. 
Their theoretical study goes back to middle nineties and starts with triangular chainsi^*^*^. 
This case attracted attention especially by the study of the emergence possibilities of fer- 
romagnetism in such systemsi^ii^, analysis of the stability of this phase^^, low temperature 
thermodynamics^, numerical studies in two band cases^, or even study of exact ground 
states in non-integrable situations^. 

The attempts to characterize the physical properties of quadrilateral chains have been 
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started later, being intensified by the study of Aharonov-Bohm cage properties^, and liave 
led to tlie description of quite interesting ground state characteristics in the diamond chain 
case which are tunable by external fields providing for example correlated half metal behavior 
applicable for spin- valve design^. 

The study of chains of pentagons intensified after 2000, when A. Heeger, A. MacDiarmid 
and H. Shirakawa obtained the Nobel prize in chemistry for the discovery and development 
of conducting polymers. The most common representatives of this class of materials are 
based on polythiophene, polypyrrole and polytriazole, all containing pentagon rings. The 
study of the physical properties of these systems has been concentrated on the search for 
flat band ferromagnetism:^ii2i>22^ and also the analysis of their ground states in exact terms 
has been started^i. 

Following the interest in increase in the number of lattice sites in the ring which forms the 
basis of the periodic chain, one arrives at chains of hexagons which are the subject of this 
paper. The hexagon chains are of real interest since are relatively abundant in nature. The 
main representative of the zig-zag hexagon chains at the level of organic materials, are the 
polyacenes. Acenes, or polyacenes are a class of polycyclic aromatic hydrocarbons made up 
of linearly fused hexagon rings. These systems attracted attention because of optoelectronic 
and electrical engineering application possibilities which have led even to organic field effect 
transistors in which pentacene is incorporated^. Furthermore, such chains present potential 
possibihties for the design of soluble acene-based transistors and circuits^, tetracene and 
pentacene have been used for the design of light emitting devices^, etc. Due to these 
properties, the polyacene structures have been studied extensively; the theoretical studies 
have been mostly confined to mean-field type of descriptions^^. One notes that the interest 
in hexagonal chains and structures is not restricted exclusively to carbon based materials, 
since also other compounds present similar cell structure, for example boron-nitride ionic 
honeycomb systems^i^ etc. 

B. The aim and technique used 

Recently, based on first principle calculations, the Coulomb interaction has been analyzed 
in acenes between two particles. It was shown^^ that the average Coulomb repulsion in these 
systems can reach values around f/g// ~ 4 — 5eV, the ratio of f/g// to the bandwidth W 



3 



being of order Ueff/W ~ 2.3 — 7.6. These results show that correlation effects are very 
prominent in this type of chains, leading the authors to the conjecture that by increasing 
the concentration of carriers in these systems, the Coulomb repulsion is expected to stabilize 
correlated magnetic ground states. 

The present paper has the aim to analyze this conjecture at a genuine many-body level. 
In order to treat the correlations properly and create premises for valuable conclusions, the 
study is made here at exact level. 

Exact results for the chains under consideration, taking into account their quasi ID 
character and non-integrable nature, are extremely rare. One only knows exact ground 
states for triangle^, and rhombus^ cases. The method used allows the exact ground state 
to be deduced under general circumstances independent on dimensionality or integrability 
and it is based on a technique which uses positive semidefinite operator characteristics. A 
positive semidefinite operator P, is the operator which has the property (0|-P|0) > for all 
vectors |0) of the Hilbert space H. It results that all eigenvalues pi of P, P\(pi) = Pi\(pi), are 
non-negative, e.g. pi > 0. Consequently, if one succeeds in finding an exact decomposition of 
the Hamiltonian H in terms of positive semidefinite terms, namely H = Pn + C, where 
Pn are positive semidefinite operators (whose number is m), and C is a constant depending 
on Hamiltonian parameters, then the operator H' = H — C has a positive semidefinite form. 
Hence, H' has a spectrum bounded below by a well defined and known number, namely 
zero. As a result, the exact ground state of H' (consequently, also of H), is given by the 
vector holding the property H'l'^g) = 0, i.e. P„|\l/g) = must be satisfied for all 
n = 1, 2, m. 

The technique which we used is based on the above presented properties. First one has 
finds an exact decomposition of the Hamiltonian of the system in a positive semidefinite 
form, and obtain the explicit form of the positive semidefinite operators P„ for all n. In 
the second step one constructs the ground state Hilbert space vector such to satisfy 
P„|\l/g) = for all n. This procedure depends on the explicit form of the P„ operators, 
but one exemplifies below a particular case which will be used for the study of the zig-zag 
hexagon chain described below. 

Let us assume that one has for n = 1,2,. ..mi < m the structure P„ = Al^An for the 
positive semidefinite operators P„, n < mi, where A„ are built up from a linear combination 
of annihilation fermionic operators. In these conditions the construction of the ground state 
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starts by the construction of the wave vector = -Bj^|0) where |0) is the bare vacuum. 
The i?^ operators (constructed from fermionic creation operators), are objects which have to 
be deduced. Their calculation is made based on the anti-commutation relation {An, -Bj^} = 0, 
which must be satisfied for all values of all indices n and fi. Indeed, if this anti-commutation 
relation holds, in -Pnl^) = ^n^ndl^ -^a')!^)' operator can be pushed in front of 

the vacuum state obtaining P„|^) = e^-^^t (H^ 5t)A„|0), which, given by the annihilation 
nature of An provides An\0) = 0. Consequently indeed -Pn|^) = holds. The phase factor 
gi<^ provides only a +1, or —1 multiplicative factor depending on the even, or odd number 
of operators in {Yl/i^ji)- 

After deducing all possible Bj^ operators, the ground state |\l'g) is obtained by restricting 
the index n E I, such to obtain \^g) = [n^ei-^illO) based on the condition P„|\E'g) = 
also for n > mi. The corresponding ground state energy becomes Eg = C. The carrier 
concentration at which |\E'g) is defined, is provided by the number of electrons introduced 
into the system by the (H^e/ -^/») operator product acting on the vacuum state |0). 

One further notes that the transformation of the starting Hamiltonian if in a positive 
semidefinite form Pn + C is not unique, can be performed in several different ways, each 
transformation places the final result in different regions of the parameter space. The pa- 
rameter space domain in which the deduced ground state is present, is fixed by the matching 
conditions. These last are relations between the H parameters and the Pn parameters which 
allow the transcription of the starting H in the used positive semidefinite form. 

The detailed presentation of the technique can be found in several recent 
publications^'^i2^'2Z. The implementation of this method at finite value of the interaction 
has been started at the end of ninetiea^^'^'^'^'^'^ and has proven to be a successful and 
powerful technique leading to exact results even in situations unexpected in the context of 
exact solutions as: three dimensions^^iSI^ disordered and interacting systems in 2D^, emer- 
gence of condensates^i^^, stripes and checkerboards in 2D^, or insulator to metal transition 
driven by the Hubbard repulsion in 2D in vicinity of half filling. We further note that 
similar techniques are used for spin models as well^. 
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C. Overview of the obtained results 

One analyzes the problem by concentrating on a fixed and given hexagon chain of zig- 
zag type (see Fig.l). The zig-zag nature must be accentuated since hexagon cells can be 
connected in a chain also in another way, namely in armchair configuration. The starting 
Hamiltonian, besides several hopping matrix elements and on-site one-particle potentials 
contains also the Coulomb repulsion as interaction, but for simplicity only at on-site level, 
providing the Hubbard interaction term with U > strength. Since we are interested in 
finding valuable information about electron correlation effects in a fixed confined system, 
one keeps the chain structure unchanged during the study, e.g. phononic contributions 
are neglected. Since Peierls transitions caused by electron-phonon interactions in chain 
structures with even number of sites per cell could influence the emerging phases only around 
half-fillingi2i'22^ one does not expect that the neglected phononic contributions will provide 
genuine changes in our results relating far from half-filling regions. 

The Hamiltonian of the starting chain is transformed first in a positive semidefinite 
form. This transformation is important not only in the context of acene structures. This 
is because several other systems of large interest today are built up from hexagons. For 
example graphene, being a 2D hexagon structure constructed from the same cell, can be 
described at the level of the transformation into a positive semidefinite form by the same 
block operators as used here for the zig-zag hexagon chain in Section IV. The block operators 
applied for the transformation were such chosen to provide a not severely restricted phase 
diagram region where the conclusions are valid. 

After this transformation one analyzes the bare band structure and one shows that flat 
bands are not possible to emerge. Hence, fiat band ferromagnetismi^ii^ (often studied in the 
context of conducting chains, see for exampl^^^^) is excluded a priori from the spectrum 
of possible magnetic phases. 

After the positive semidefinite form of the Hamiltonian has been obtained, one turns 
to the construction of the ground states l^^g). This step starts with the deduction of the 
operators denoted hereafter by Bj^, building up l^^g)- One notes that the label index has 
also a spin component. Interestingly, it turns out that all these operators are extended, i.e. 
have components placed in each cell and cannot be introduced in a restricted domain of the 
chain with extension smaller than the system size in the direction of the primitive vector. 
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The extended nature of all IBj^ operators is a special characteristic of the hexagon chains, 
which has not been observed previously in the case of triangle, quadrilateral and pentagon 
chains^i^i^i. This aspect is important for two reasons. First, the methods for the treatment 
and deduction of the extended operators in the context of the positive semidefinite operator 
technique are practically completely open and unknown. On this line one notes that even 
if it has been shown that the treatment in another context of the extended operators is 
feasible^, one has only one case described quite recently in the literature, namely in the 
study of the insulator to metal transition driven by the Hubbard repulsion in 2D in the 
vicinity of half filling^^. Since a strategy developed for 2D, directly cannot be applied for 
quasi ID structures, at this point a special technique has been constructed and described 
for the deduction of the Bj^ operators, which can be used for other chain structures as well. 
Second, the extended nature of Bj^ terms is present in a confined space region delimited by 
the chain itself. Consequently, by increasing the carrier concentration (i.e. increasing the 
number of independent Bj^^ operators in l^^g)), the Bj^^ operators present in the ground state 
wave vector will satisfy connectivity conditions (i.e. will act on common lattice sites). Given 
by this, the spin indices of the electrons must be correlated in order to avoid as much as 
possible the double occupancy and to reduce in this manner the ground state energy. This 
is the route for the emergence of ferromagnetism in these systems. The obtained results 
show that the conjecture by Brocks et al.— made at two-particle level for acenes is indeed 
viable in a rigorous many-body frame provided by a fixed zig-zag hexagon chain described 
by a Hubbard type of model. 

The remaining part of the paper is structured as follows: Section II. describes in detail 
the system of interest and its Hamiltonian, Section HI. analyzes the non-interacting band 
structure. Section IV. presents the transformation of the Hamiltonian in a positive semidef- 
inite form. Section V. solves the matching conditions. Section VI. presents the deduction 
strategy for the exact ground states, and the obtained ground states together with their 
physical properties, and Section VII. containing the summary and conclusions closes the 
presentation. The Appendix contains mathematical details related to the exemplification in 
a simple case of the deduced ground states. 
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II. THE SYSTEM OF INTEREST 



The zig-zag hexagon chain is presented in Fig.l. In the cell constructed at the site 
i one has 4 lattice sites whose position relative to the site i is given by the vectors r^, 
z/ = 1, 2, 3, 4. The system contains four sublattices S^, containing the lattice sites i + r^ E S^, 
hence u represents as well the sublattice index. For mathematical simplicity one considers 
ri = during the calculations. The primitive Bravais vector a is directed along the line 
of the chain. The hexagon side is b = a/-\/3, where a = |a| is the lattice constant. One 
considers below a chain build up from Nc cells, the number of lattice sites being A^a = 4:Nc, 
while the number of electrons is denoted by A^. 

The Hamiltonian of the system H = Hq + Hjj is given by 

a i=l 

+ i^^lc[_|_ri,(7^i+r3,o- + H.C.] + \t (Ci4_r2,cr^i+i"3>o' ^i+r3+a,(T^i+r2,o- + c|_^r4,(T^i+i"i+a,(j 

^i+ri,cr^i+r4,o-) + ^1 (Ci^ri+a,o-^i+ri,cr + c|_|_r3+a,o-^i+r3,CT) 

+ '^e(c|_,_j.2+a,(T^i+r2,o- + c|_,„r4+a,(T^i+r4,cr) + H.C.] 

+ [fo(^i+r2,o- + T^\+VA,a) + Cl(^i+ri,(T + '^i+rs.o-)]) 
Nc 4 

^t/ = J]]^^i+r,.,T^i+r.a) (1) 
i=l y=l 

where cj^ is the canonical Fermi operator creating an electron at the site j with spin projec- 
tion 0", and fii^a = cl^Ci^cr represents the particle number operator for electrons with spin a 
placed at the lattice site i. Since real systems are in view, several hopping matrix elements 
are considered in ([T]). The t and ti terms characterize nearest neighbor hopping matrix 
elements {ti is placed perpendicular to the line of the chain). The t', t[, contributions 
represent next nearest neighbor hopping matrix elements {t'l and are directed along the 
line of the chain, t'^ connects internal sites while tf. external sites of the chain). The eo, 
ei parameters are on-site one-particle potentials (eo is placed on external sites), see Fig.2. 
Finally, Hu characterizes the on-site Coulomb repulsion where ?7 > is considered. 
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FIG. 1: The studied zig-zag hexagon chain. The ceU constructed at the site i contains 4 sites whose 
in-cell positions relative to i are specified by the vectors r^,, v = 1,2,3,4, where, for convenience 
one takes ri = 0. The vector a represents the unique primitive vector of the Bravais lattice. In the 
same time, denotes lattice sites in four different sublattices 3,^. 




FIG. 2: The parameters of the kinetic part of the Hamiltonian. The t and ti terms represent 
nearest neighbor hopping matrix elements from which ti is placed along the touching bonds between 
hexagons, t', t'l and te are hopping matrix elements describing next nearest neighbor hoppings. 
From these te describes the external hopping relative to hexagons, while t[ is parallel to the line 
of the chain, being placed inside the hexagons. Finally, ei, (eo) represents the on-site potential on 
contact points between hexagons (external sites of hexagons). For the clarity of the notations, the 
next nearest neighbor hoppings are separately presented in the third cell of the figure. 
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III. THE NON-INTERACTING BAND STRUCTURE 



The deduction of the bare band structure 



In order to deduce the non-interacting band structure one transforms from ([T]) in 
k-space. For this one Fourier transforms the cj.o- operators as follows 



^^'^ k=i 



-ik(i+r^+r) 



(2) 



where Cj^^ico- represents the annihilation operator for the state (k, a) in the sublattice u. Note 
that one has ri = 0, and r from ([2]) takes two possible values, namely and a. Substituting 
(E]) in (H]), one finds 

o- k 

+ [i(:i,4(k)c[,^^^C4,k,a + t2,3(k)4,k,<7^3,k,<7 + ti,i(k)cJ,k,aCi,k,a 

+ t2,2(k)c^ ,^C2,k,<7 + t3,3(k)4,k,aC3,k,a + t4,4(k)c^_,^ ,^C4,k,<7 + H.C.] 

+ [eo(4,k,aC2,k,a + 4,k,aC4,k,<x) + Cl (4 ,k,aCl,k,a + 4,k,aC3,k,a)] } , (3) 

where the first, second and third, and fourth rows represent in order the nearest neighbor, 
next nearest neighbor, and on-site contributions. Furthermore one has 

ti,2(k) = t(e+^'^("-'"^) + e"''^'"^), ti,3(k) = t.e-'^''^, ts^k) = t(e+''^('"^-'"*) + e-^^^^''^--^-'^)), 

ti,i(k)=t;e+^^ t2,2(k) = tee^^'^, t3,3(k) = t;e+^'«^, t4,4(k) = t^e^^^ (4) 

where for r^, v = 2,3,4, see Fig.l, while a represents the primitive vector of the Bravais 
lattice. Taking into consideration a Cartesian system of coordinates with versors (ii,ji) 
whose X-axis is directed along the line of the chain, one has 



a a . a , a. 

where a = |a| holds. 

Introducing the 1x4 row vector and its 4 x 1 adjoint C by 

C2,k,cr 
C3,k,(T 
V C4,k,,7 / 



a\/3 . 



(5) 



~ (^l,k,o-' '^2,k,o-' '^3,k,cr' '^4,k,(T;; 



(6) 
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one observes that ([3]) can be written in a matrix product form, namely 



(7) 



cr k 



where, introducing the notation /(k) = (1 + e"*"**^), for the matrix M one has 



M 



/ei + 2t;cosak te-*'^'"V(k) t^e-'^'''' t'e~'^''* f{k) \ 

te+i^^^f*(k) eo + 2tecosak t'e-*^(''»-''2)/*(k) 



(8) 



^^g+ikra ^/g+ik(r3-r2)j(^]^) ei + 2t'i COS ak te-*''(''4-'"«)/(k) 
\ t'e+'^''^f*{k) te+*''(''4-'-3)/*(k) eo + 2teCosak / 



The bare band structure is obtained from the secular equation of the matrix M from 
which provides 



(eo - A + 2te cos ak)2 [(ei - A + 2t[ cos ak)^ - tj] + (eo - A + 2t^ cos ak) |/(k) p 

x[2ttit' - (ei - A + 2t'iCosak)(t'^ + - t\f{'k)\^[t{eo - A + 2te cosak)(ei - A + 2ticosak) 

-(eo - A + 2teCosak)t'ti + \f{k)\H{t'^-t^)]~t'\f{k)\'^[t'{eo - A + 2teCosak) 



The non-interacting band structure containing 4 bands is obtained from the solutions 
Er,{k) = A of the quadratic algebraic equation ([9]), which provides four solutions r] = 1,2, 3, 4. 
One further notes that in |/(k)p = 2(l + cosak) holds, and = ak G (— vr, tt] is satisfied 
for the first Brillouin zone. 

Introducing the notation ia = ea — A for a = 0, 1, from ([9]) one obtains the equation for 
the bare band structure in the form 



B. The properties of the bare band structure 

A picture representing an exemplification for the bare band structure is presented in 
Fig. 3. Modifying the parameter values, one obtains cusp points which can be considered in 
a given extent reminiscent of the Dirac points present in the graphene case^, see Fig.4. But 



x(ei - A + 2t;cosak) - (cq - A + 2te cos ak)tti + \f(k)\H'{t'^-t'^)] = 0. 



(9) 



(eo + 2teCosA;)^[(ei + 2t[coskf - tl] - 4(eo + 2teCosA;)(l + cos A;) 
X [(ei + 2t[ cos k){t'^ + t'^) - 2ttit'] + 4(1 + cos kf{t'^ - t^f = 0. 



(10) 
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FIG. 3: An exemplifying image for the band structure plotted on a 27r domain for the k variable. 
The used parameters are h/t = l,t'/t = 0.2,t[/t = 0.2, te/t = 0.2,eo/t = l,ei/t = 1. By 
modifying the Hq parameter values and signs, the relative inter-band distances and band shapes 
can be changed. 




ak 

FIG. 4: Band structure with cusp points plotted at ti/t = 1.8/1.6, t' /t = 0.5/1.6, t'^/t = 
0.6/1.6, te/t = 0.2/1.6, eo/t = 1, e\/t = 2/1.6 parameter values. 

since in the studied case, the \k ± k*\ values not provide the same energy around the cusp 
points k*, for small k, Dirac points in rigorous terms are not present here. Another case 
with closely situated intersection points for the lower two bands is presented in Fig. 5. 
It is important to underline that flat bands are not possible to occur for the studied 
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FIG. 5: Band structure with closely situated intersection points plotted at ti/t = 1.4, t'/t = 
Q.2,t'i/t = 0.2,te/t = 0.2,eo/t = 3, ei/t = 7.76 parameter values. 

system. In order to show this one writes (fTO|) into the form 

A4 cos^ k + As cos^ k + A2 cos^ k + Aicosk + Aq = 0, (11) 

where An = An{eo,ei,t,t' ,ti,t[,t(,, X) for all n = 0,1,..., 4 are k independent. Flat bands 
(e.g. k independences in A) are obtained if simultaneously An = holds for all n. This 
provides the following system of equations 

Asieo, ei, t, h, t', t[, t^, A) = 16(eo - A)tetf + 16(ei - \)tlt[ - I6tet[{t^ + t'^) 

A2(eo, ei, t, ti, t', t[, t,, A) = 4tf (eo - A)^ + 16(eo - A)(ei - X%t[ + AtHe^ - - At^tl 

-met[{t^ + t'^) - 8t[{eo - A)(t^ + - 8te(ei - A)(t^ + t'^) + Wttit'te + 4:{t'^ - t^f 
Ai(eo, ei, t, ti, t', t'l, te, A) = 4(eo - A)'(ei - X)t\ + 4(eo - A)(ei - Xft^ - At%{e^ - X) 

-8t;(eo - X)(f + - 8te(ei - X){e + t'^) - 4(eo - A)(ei - X){e + t'^) 

+8ttit'(eo - A) + mtit'te + 8(t'^ - t^ f 
Ao(eo, ei, t, ti, t', t;, te, A) = (eo - Xf^e^ - Xf - tl^e^ - Xf - 4(eo - A)(ei - X){e + t'^) 

+8ttit'{eo - A) + 4(t'2 - ef. (12) 

As seen, (fT2l) presents An = solutions for all n only at tet'i = 0, hence taking non-zero 
values for all introduced hopping terms, fiat bands are not possible to occur into the system. 
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FIG. 6: The lattice sites on which the block operators A^^i ^j, u = 1,2,3,4 are defined for the cell 
constructed at the site i. In each figure the position of the site i is explicitly shown. The numbers 
are representing the numbering of the coefficients of A^,i,a operators (see also ()13p ). 



IV. THE TRANSFORMATION OF THE HAMILTONIAN IN POSITIVE 
SEMIDEFINITE FORM 

One introduces four block operators Ai^^i^^r, v = 1, 2, 3, 4 defined in the cell constructed at 
the site i by the following relations 

^l,i,(T = fllCi+ri,(T + 02Ci+r2,o- + C^sCi+rs.cr + '3oCi_a+r2,o-) 

^2,i,o- = &lCi+ri,(T + &3Ci+r3,CT + &4Ci+r4,o- + &0Ci-a+r4,o-, 

^3,i,(T = 'J^2Ci+r2,(T + C^sCi+rs, cr + C^oCi+a+rs.cD 

^4,i,(T = ClCi+ri,o- + e4Ci_(.r4,o- + CoCi+a+ri,o-) (13) 

The lattice sites on which the operators Ai/,i,<T are defined for the cell constructed at the site 
i are presented in Fig. 6. 

Using periodic boundary conditions and based on ( fT3|l . the Hamiltonian of the problem 
transforms into a positive semidefinite form as follows 



Nc 



(14) 



i=l a u=l 

where N denotes the operator of the total number of particles. The relation (fH|) holds if 
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the following matching conditions are satisfied 



t = 


* * 

— 02(31 — Qj^CLq — 


= &4^3 = blbo, 






te 


= 02(30 = b^bo, 








tl 


= O3O1 + &3&1, 








t' 


= O2O3 + d^d^, = 


= O3O0 + (io(^2 


= 6461 + 6461 = 


= blbo + 6064, 


t[ 










eo 


- is: = loop + 


|a2r + = 


\bo\' + M' + 






-K= |oip + 


|eor + |eip + 


\bi\' = \b3\' + 





(15) 

The matching conditions f|T5l) have been obtained by calculating J2a S^=i ^lia^uXa 
from f|T^ and equating the result term by term to the expression of Hq from ([1]). 

V. DEDUCTION OF i^i <^ OPERATORS 

In deducing the block operators defined in f|T3|) one must calculate the prefactors 
aa,ba,da,ea prcscut in equations (fT3l) . These are obtained by solving the matching con- 
ditions ([T5|) . 

Solving the equations present in (|T5l) one proceeds as follows. The first two rows of Eq.( JT5l) 

provide solutions only for te > and taking into account real hopping matrix elements, give 

t 



|oi = \b: 
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The fifth row of (fT5|) provides 



a2 = ao = — , 64 = Oo = 7T- (16) 



eo = ^, c/o = ^, (17) 



and the remaining equations from f|T5l) transform into the form 



Oi Oi CX3 

./ ,^1 , * bi el 
t =t— + 6461 = t— + — , 
&3 &3 el 

eo - if = 2te + = 2te + |e4p, 

.2 ,2 fi2 

ei - ^ = ^ + + |eiP + l&iP = ^ + ^ + |rf3p + |a3p. (18) 

te |ei|^ te Idgr 
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The second and third hne of (fTSi) arises from the fourth row of (fTSi) . For the first hne of 
f|T8|) . given by |ai| = I63I in f|T6l) . one uses 63 = aie*'^?. Now from the second and third hne 
of (ITS]) , besides the condition t'l > 0, one finds {d^l = |ei| = \/t[. Consequently, the last 
two equations from flTHl) give \d2\ = {e^l, and = [03 1. Using bi = 036**^3 from the previous 
condition, the first equation of (ITS]) becomes 

ti = a*ai + a*a3e'('^3-0?), (I9) 

which, because of the real nature of ti gives ^3 = (pl, from where ti = a^ai + 0*03 = 
2i?e(a*a3), furthermore 61/63 = 03/01 arises. At this step one observes that in (1131) each A^ i ^^ 
operator can be multiplied by an arbitrary phase factor without changing the expression of 
the transformed Hamiltonian (tH]). Consequently, by this multiplication relating Ai^i^cr, the 
prefactor oi can be taken real, hence from the expression of ti the parameter 03 must also 
be real. As a consequence, ti = 2aia3 holds, from where 

As a result, one finds that the coefficients are given by 

/— t I— , , 

OO = V^e, «1 = —?=■, «2 = V^e, «3 = cy, ■ (21) 
V 

In what follows one concentrates on the coefficients. Given by h\ = 036**^3^ 63 = aie**^! 
where ai, are real and 0" = ^3, a multiplication by e^**^" of v42,i,(j from (IT^ provides real 
bi and 63. Hence ( 1T6H2T]) together with 61/63 = a^/ai = tite/{2t^) gives 

60 = v^, &i = ^, ^^ = ^- (22) 

At this moment one continues the calculation by deducing Cq, and da- From the real 
value of if:', taking into account the second and third expression of (ITS]) , one observes that 
by fixing ^3 (similarly ci) to be real by multiplication by a phase factor, then ^2 (similarly 
64) also becomes real. Consequently t' = (tite)/(2t) + d2\/tf[ = (tite)/(2t) + 64 a/^^ and 
dz = ei = A/t^ together with Eq. (IT7]) provides 

2tt' — i\te 



eo = V^'i, ei = Vt'i, 64 



rfo = v^, ^2 = ^^^^^, rf3=A. (23) 
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From the last two lines of (fTSi) one further finds the K value 

{2U' - ^ ^ 

K = e,-2U- ^ ^^^^^ , (24) 

and a condition which must be satisfied by the Hamiltonian parameters, namely 

eo = e, + 2(t.-t;)-^-§ + i^^^^^. (25) 

The equality fl25l) represents the requirement for the on-site potential eo necessary to be 
satisfied in order to find the expression for the transformed Hamiltonian (IT^ a valid relation. 
Besides fl25]) . the parameter space region where the transformation from ([T]) to the positive 
semidefinite form (JT^ of the Hamiltonian is valid, is described by > 0,t[ > 0. In the 
presence of all these conditions, the bare band structure is exemplified by Fig.5. 

Using the results presented in fl2m23p one finds for the block operators introduced in f|T3|) 
the explicit expressions 



^l,i,(T V ''el^Ci+r2— a,(T T~ ^K-'rVz^cy) ~r at— Ci-fri,(T ~r Cj+rsjcrj 



^2,i,CT V ^e(^Ci+r4— a,(T ~r Ci+r4,(Tj ~r ^^y^Ci-|_r3^o- ~r 2^ Ci+ri,(T5 
1 FT / " " \ 2tt' — t-\t^ „ 



VI. THE DEDUCTION OF EXACT GROUND STATES 

The ground states are deduced in two steps, i) First one looks for 5^ operators proper 
for the construction of the ground state by 

i*.)=n^»' (27) 

where |0) is the bare vacuum, / is a set of parameters /i, and the Bjj^ operators satisfy 

{i,,i,,,5t} = 0, (28) 

for all values of all indices z/, fi, a and i. The motivation for this step is that given by the 
requirement (J28l) . for the wave vector |^) = n,.i -^^lO)' '^'^^ property 

EEE^UAi,.i^) = o- (29) 



i=l a v=l 
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FIG. 7: The coefficients Xi^^i from (j30p of the sites of the hexagon chain presented in the neighbor- 
hood of the lattice site j. 

Indeed, ( 129|) holds since starting from fl28|) . the A^^i^„ operators in (129!) can be pushed in 
front of the vacuum state, hence given by A,^_i o.|0) = 0, the equahty in Eq. fl29l) is satisfied. 

In the second step ii) one chooses from all the possible Bj^ operators those Bj^^, I = 
{/ii,/i2, .••,/iAf}, which satisfy also ^^i/dl^^e/ ^/llO)) = 0- 

A. The operators needed for the construction of the ground state 

1. Requirements for the operators building up the ground state 

One considers the -B^ operators as the most general linear combination of creation op- 
erators with fixed spin projection acting on each lattice site of the system. Consequently, 
besides a fixed /i = /ii index, also the spin projection a will be explicitly used as a label, 
and one has 

j 

where the site j runs over all lattice sites of the sublattice 5*1, and the x^^i coefficients 
represent the unknown variables which must be deduced. The placement of the coefficients 
x^^i is presented in Fig. 7. 

Now using (I26f30l) in one finds the following linear equations for the x^^i coefficients 



18 



from (l30i) defined for the cell constructed at the site i = j 



(31) 



where the notation Q = (2tt' — tite)/ {2t^t']) has been introduced. From (!3T!) one finds 



Qt Q 

f^3,i H 7=2;4,j-i, 



6p \ / 6 



X2, 



t ti 



(32) 



The equations from (152]) can be written in a matrix form as 

where is a column vector with 4 components, which taken in order are 

given by xij+i, xsj+i, X2j, a;4j. Similarly, Zj_.i is the column vector with elements 



Z,_i 



Finally, the 4x4 matrix R is defined by 



a;3,j 
V ^4,i~l / 



(34) 



R 



- 1 







Qt Qtl _ ^ _Q_ g 

t *i _i 



*1 

2t 



v 



*1 

2t 







(35) 



Starting from (153]) . in the following subsection one shows how it is possible to construct one 
(starting) -B^^ o- operator. 
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2. The construction of one starting Bji^^„ operator 
Based on (l33i) one finds 

Zj+m = W^Z,, = (R)'". (36) 

Consequently, since periodic boundary conditions are taken into account (i.e. after a finite 
number of m steps, in Zj+m one must recovers Zj), an operator B'^^^^ from fl30|) . proper 
for the construction of the ground state, is obtained when the matrix has at least one 
eigenvalue 1. In this case the Xa^i unknown parameters from ( l30i) are given exactly by the 
corresponding Z| eigenvector 

W^Z] = Z^^. (37) 
One finds in this manner the starting coefficients Xa^i in the form 

For the next (m — 1) cells Z^/, the coefficients xij/, xsj/, X2j'-i, a;4j/_i are obtained by suc- 
cessive application of the R matrix to Z^ as follows: Xij+i, Xaj+i, X2j, a;4j are obtained 
from Zi deduced as Zi = RZ|; a;ij+25 2:3^+2, a;2j+i, a;4j+i are obtained from Z2 deduced as 
Z2 = (R)^Z|; etc.; Xij+^-i, a^sj+m-i, a;2j+m-2, a;4,j+m-2 are obtained from Z^^i deduced 
as Zra-i = (R)™"^Z|. Given by ([37]), for the coefficients Xij+rn,X3j+m,X2,j+m~i,X4j+rn-i 
entering in Z^ = (R)"Z^, the values xij+rn = xij,X3j+rn = xsj, X2,j+m-i = 
X2,j-i, X4j+m-i = X4j_i are reobtained. After this step, starting from the coefficients 
xij+m+i, X3j+m+i, X2,j+m, X4j+m prcscut in Z^+i, the expressions are periodically repeated 
since Z^+i = R-Z^ holds again. In this manner a Bj^_^ ^ operator has been constructed for a 
system characterized hj = p x m, where p is an arbitrary positive integer. 

3. Generation of new bJi^ ,0- operators 

Once a 5^^^ ^ operator has been deduced following the procedure described above, one 
has a proper operator for the construction of the ground state that has the obtained Xa^i 
coefficients structure holding a periodicity of m Z cells, the starting cell being the jth 
cell Zj. If one moves the starting cell to the Z^+i cell, one obtains a new, (usually) linearly 



20 



independent ^ operator. Repeating the procedure, starting from ^, maximum (m— 1) 
new operators Bj^^^^, (3 = 2,3,4, ...,m can be created (the hnear independence must be 
checked at each step). If more than one unity eigenvalues are present, the procedure can 
be repeated for each eigenvalue, constructing in this manner for the fixed m periodicity, 
maximum qxm different ^ operators, where q represents the number of unity eigenvalues 
corresponding to eigenvectors with non-zero norm. An exemplification is presented in the 
Appendix. The procedure can then be repeated for another m. 

4- The ground state wave function 

Using the operators deduced above, the |\E') wave vector described below (!28|) has the 
form 

i^)=nn^uio)- (39) 

M Cm 

One notes that the Bj^ ^ operators (given by fl33l) ) are all extended operators. Hence increas- 
ing the number of carriers and taking into account that all 5^ ^ operators are confined in 
the same quasi ID structure, it results that connectivity conditions will be satisfied between 
different Bj^^ ^ operators with different /i indices (e.g. Bj^^. ^ and Bj^^ ^ operators for i ^ j 
will touch each other, hence will act on a given finite number of common sites). For the 
example provided in the Appendix A. this happens starting from the N = 2 number of 
particles. Consequently, in order to minimize the ground state energy, all a indices in ( 139|) 
will be fixed to the same value (a^ = 0"^/ for all fi, fi') providing in this manner zero double 
occupancy, hence a ferromagnetic ground state of the form 

i^.)=n^icio)' (40) 

where a is fixed. Indeed from (140|) . as constructed from (!39|l . will satisfy 

iJ^iti J2a Ylt=i ^li o-^i'.i.o-)!^^) = for the first part of the transformed Hamiltonian from 
( Tn|l . while Hul'^g) = will be given by the absence of the double occupancy. Consequently, 
|\Efg) represents the ground state of H from (HM . the ground state energy being Eg = KN, 
where K is given in (12^ . Since one Bj^^^ operator introduces in fact one particle into the 
system, the particle number (hence the concentration of carriers) at which |\E'g) in (1401) is 
defined, is given by the number of components of the product from (HOj) . 
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As described, represents a ferromagnetic state. This ground state is itinerant and 
metallic because of the following reasons: i) the ground state wave function is extended, and 
ii) up to the N^ax value equal to the number of components of the product from (HUi) . one 
has for the particle number (A^ < N^ax) dependent chemical potential /i^, the expression 
Sfic = /ic(+)-/U,(-) = where = Eg{N+l)-Eg{N) and/i,(-) = Eg{N)-Eg{N -1), 

see also Ref.f^]. 

Our up to date results show that the ferromagnetism emerges in the low concentration 
domain, a pedagogical example being presented in Appendix A. If one characterizes the 
carrier concentration hj n = N/Nc (i.e. average electron number per cell), the exemplified 
case in Appendix A. of six explicitly given operators describing through (1401) six interacting 
electrons, even for a chain made of 36 cells (for example) gives n = 1/6, a far from zero finite 
concentration value. Consequently, the reported results clearly demonstrate the presence of 
itinerant ferromagnetism in the systems under study at low concentration, for finite chains 
(treated with periodic boundary conditions). Since 36 cells is a huge number for oligo- 
acenes (i.e. few linearly fused hexagons) placed in the center of the attention^^, our findings 
represent genuine information for potential application possibilities. 

Furthermore, one underlines that since (as shown in Sec. III.) flat bands are not present 
in the system, the here described ferromagnetism is not of flat band type, and is provided 
by a joint effect of correlation and conflnement. 

One notes that for the flnite sample case explicitly exemplifled in Appendix A, the Hub- 
bard repulsion must satisfy only the U > condition, the U value itself being without 
importance. The phase diagram region where the presented solution occurs is not severely re- 
stricted. Besides two sign requirements present for two hopping matrix elements {te, t'l > 0), 
only one condition for on-site one-particle potentials is present (see fl25l) ). which requires 
a flxed value for eo acting on external sites. Since this on-site potential can be tuned by 
an external gate potential^., the described itinerant ferromagnetism can be in principle even 
switched on by an external electric fleld. 

A detailed study of the phase diagram and the behavior in the thermodynamic limit is 
in preparation, and will be published elsewhere. 
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VII. SUMMARY AND CONCLUSIONS 



A zig-zag hexagon chain described by a Hubbard type of model containing on-site 
Coulomb repulsion is analyzed in exact terms by a technique based on positive semidef- 
inite operators. The calculations effectuated with periodic boundary conditions provide 
exact ground states of ferromagnetic and itinerant nature in the low concentration region. 
Flat bands are not present in the non-interacting band structure, and the ferromagnetism is 
created by a joint effect of correlations and confinement. The parameter space region where 
the described phase emerges is not severely restricted, and the unique on-site one-particle 
potential which must have a fixed value can in principle be tuned by external gate potentials. 

The reason for the emergence of ferromagnetism is connected to two aspects. First, all 
operators entering in the construction of the ground state turn out to be extended. This 
is an interesting property of the studied system since it has not been observed for triangle, 
quadrilateral, and pentagon chains^i^i^i. Second, all these extended operators with contri- 
butions in each cell act on lattice sites confined in the chain structure under consideration. 
Hence by increasing the carrier concentration, connectivity conditions between the operators 
of the ground state wave vector will emerge. These connectivity conditions (i.e. different op- 
erators describing different electrons act on common lattice sites) fix the spin indices to the 
same value, minimizing in this manner as much as possible the double occupancy in order to 
reduce the ground state energy, and provide the ferromagnetic ground state. The itinerant 
nature of the obtained ferromagnetic ground state may open new application possibilities in 
spintronics. 

One adds below some observations related to the exact nature of the results. One knows 
that the flat band, or nearly fiat band ferromagnetism can be described in exact termsF^^. 
Furthermore, it was shown that ferromagnetism emerging for completely dispersive band 
structure provided by spin-spin interactions in itinerant systems can also be described in 
exact terms^. The here reported results enlarge the spectrum of exact descriptions relating 
ferromagnetism in completely dispersive systems: one demonstrates that when the joint 
action of confinement and correlations is present, itinerant ferromagnetism can be described 
in exact terms at finite value of the interaction even in the frame of a non-extended Hubbard 
type of model. 

The deduced results deserve a further observation as well: In the search for organic 
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ferromagnetism is no more necessary to look for hopping and on-site one-particle potential 
values in different chain structures in order to introduce flat bands in the system leading 
to flat band ferromagnetism. Itinerant ferromagnetism can be obtained even in completely 
dispersive systems, for example by a joint effect of correlations and confinement. 



Acknowledgements 

Z.G. gratefully acknowledges financial support provided by the Hungarian Research Fund 
through Contracts No. OTKA-T48782 and the Alexander von Humboldt Foundation. 



APPENDIX A: EXEMPLIFICATION FOR THE CONSTRUCTION OF B^^ OP- 
ERATORS 



Let us consider the following pedagogical example for the construction of operators needed 
in the ground state wave vector. Using the notations a = Q/ ^/t[, 6 = t/tg, c = ti/ (2t), the 
R matrix from fl35l) becomes 



/ an — 



R 



ac — 1 ah a \ 
ah ac — 1 a 
-h -c -1 
-c -h -1 / 



(Al) 



Using for simplicity the conditions c = —6, a = —1/6^, for W4 one obtains four unity 
eigenvalues. One of eigenvectors has only zero elements (i.e. zero norm), so this cannot be 
used to construct physical states, but three other eigenvectors have non-zero norm. These 
are the following ones 

/o\ 

1 



/ \ 



-1 
V 1 / 












(A2) 
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Let us consider first the eigenvector. For it one has 

/ \ 





RZ^ 



1 



/ i ^ 



R^Z^ 



-1 
V 1 / 



R^Z^ 



/ \ 



-1 

V 1 y 



(A3) 



Taking Zj = Z^ one finds based on ( ]A3I) the following starting -B^^ q. operator 

••• + (-4+r„. + cJ+r,J + (-Tc| 



If If t t 

— r' -I r' -\- r' — r' 

+ (^j+2a+r2,cr ~ '^j+2a+r4,cr) + (■^'^j+4a+ri,cr ~ ■^'^j+4a+r3,o- ~ '^j+3a+r2,o- + '^j+3a+r4,o-) 



J 



+ -c 



+ C: 



(A4) 



'"j+4a+r2, cr ^ '"j+4a+r4,o 

where the contributions from fixed Z cells are collected in parentheses. The second fi^^ ^ 
operator is obtained from (]A4|) by translating the starting point of the period with one cell, 
i.e. taking Zj+i = Z^ One gets 



...+ tC- 



— c 



+ c 



t 



+ -c 



+ c 



t 



^2,0- ••• ' \^^j+a+ri,CT ^^j+a+ra.o- ^}+r2,a ' "^j+r4,CT/ ' V ^j+a+r2,CT ' '^j+a+r4,o- 



1 



^■^j+Sa+ri.cr ^ ^ "-j+Sa+rs.o- ^ '-j+2a+r2, cr "-j+2a+r4,cr; V^j+3a+r2,o- ^j+3a+r4,(7; 

(A5) 



(^4+5a+ri,cr ^ '^j+Sa+rs.cr 4+4a+r2,o- + 4+4a+r4,o-) + • • • • 



A further translation of Z^ to the following cell by taking Zj+2 = Z^ provides — -B^^ o- while 
the fourth translation by taking Zj+a = Z^ reproduces i?^^ ^. Hence these steps do not 
provide new linearly independent ^ operators. 

Now one concentrates on the Z| eigenvector. In this case one obtains 



RZo 





b 

\-bJ 








R^z:^ 



-1 
-h 

\ b J 



R^Z^ 



1 


voy 



(A6) 



Taking Zj = Z2 one obtains the third linearly independent ^ operator as follows 



Bl. „ = ... + (4,,+r3,J + i-^U.^- + bc\ 



2a+ri,o- ' "^j+a+r2,o- "^j+a+r4,CTV 



(^j+3a+ri,<T) + ( ^jH 

,t 



■4a+r3,(7 



j+3a+r2,(T 



j+3a+r4,(T/' 



-j+5a+r3 



(A7) 
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The fourth i?^^ ^ operator is obtained from flA7p by translating the starting point of the 
period with one cell, i.e. taking Zj+i = Z|. One gets 



+ 
+ 



-C 



t 

j+3a+ri,o- 

t 



^+2a+r4,<x) + 



''^j+5a+r3,(7 ^^j+4a+r2,cr + ^^j+4a+r4,cr) + (^j+ea+rs.cr) + 



-j+4a+ri,(T/ 
t ^ 



(A8) 



Continuing the translation by one cell and taking Zj+2 = Z2 one obtains the Bjj^^ ^ operator 
as follows 



S' = + fr' ) -\- f-r' - ftr"'' + he' ) 



^ l'-j+3a+r3,o-/ ^ V ^j+4a+ri,o- ^ "^j+3a+r2,o- ^^j+3a+r4,cr J 



(A9) 



With Zj_|_3 = Z2 one obtains the last independent operator generated by Z|, namely -B^g^o-, 
in the form 



Mere- 



st 



+ ( ^j+a+ri,o- + ^^j+r2,o- ^^j+r4,o-) + (^j+2a+ri ,o-/' 



+ 
+ 



-c 



t 

j+3a+r3,(T 



^Cj+2a+r2,(T + ^^j+2a+r4 



J + (4 



"'^j+5a+ri,(T + ^4+4a+r2,cr 



^4+4a+r4,o-) + (^j+6a+ri,cr) + 



-j+4a+r3,(T; 
t N 



(AlO) 



The use of Z3 leads to 5^^^ operators linearly dependent on -6^3,0-5 -^^4,0-' -^Ate.o-- 

Further operators 5^ from the eigenvectors of W4 presented in (]A2p cannot be con- 
structed. But the list of new independent Bj^^ operators can be continued by looking for 
unity eigenvalues of W„, n 7^ 4, and continuing the procedure presented above. For exam- 
ple in the present case W2 has also unity eigenvalue holding eigenvector with non-zero norm 
Z| with properties 

-1 



1 




RZ4 





V / 



(All) 



leading to following contribution by Z, 



26 



The linear independence of the exphcit ]3j^ ^ operators present in (lA4|A5IIA7l - IA10p results 
from the fact that each of them has a specific (and different) group of missing initial fermionic 
creation operators. 
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